Direct measurement of quantum phases in graphene via photoemission spectroscopy 
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Quantum phases provide us with important information for understanding the fundamental prop- 
erties of a system. However, the observation of quantum phases, such as Berry's phase and the 
sign of the matrix element of the Hamiltonian between two non-equivalent localized orbitals in a 
tight-binding formalism, has been challenged by the presence of other factors, e. g. , dynamic phases 
and spin/valley degeneracy, and the absence of methodology. Here, we report a new way to di- 
rectly access these quantum phases, through polarization-dependent angle-resolved photoemission 
spectroscopy (ARPES), using graphene as a prototypical two-dimensional material. We show that 
the momentum- and polarization-dependent spectral intensity provides direct measurements of (i) 
the phase of the band wavefunction and (ii) the sign of matrix elements for non-equivalent orbitals. 
Upon rotating light polarization by n/2, we found that graphene with a Berry's phase of mv (n = 1 
for single- and n = 2 for double-layer graphene for Bloch wavefunction in the commonly used form) 
exhibits the rotation of ARPES intensity by ir/n, and that ARPES signals reveal the signs of the 
matrix elements in both single- and double-layer graphene. The method provides a new technique 
to directly extract fundamental quantum electronic information on a variety of materials. 

PACS numbers: 79.60.Jv, 03.65.Vf, 31.15.aq, 81.05. uc, 73.22. Pr 



I. INTRODUCTION 

Quantum phases are the most beautiful example of 
quantum physics and essential to understand physics in 
any material. For example, Berry's phase, the accumu- 
lated phase in the eigenfunction acquired by evolving the 
quantum system adiabatically around a closed loop in the 
parameter space of the Hamiltonian [1], has been shown 
to be responsible for the Aharonov-Bohm effect [2J, the 
half-integer quantum Hall effect [MS], etc. Another im- 
portant example is the sign of the hopping matrix ele- 
ment (or hopping integral) (<fii \ H \<p2) of the Hamiltonian 
between two non-equivalent localized orbitals (j>\ and 4>2 
in a tight-binding formalism. This phase, a fundamen- 
tal quantity in determining the electronic structure of a 
system, is dictated by the characteristics of the atomistic 
interaction, e. g. , whether it is attractive or repulsive. 
Both of these are important to directly extract funda- 
mental quantum electronic information on a variety of 
materials [SHE]. 

In graphene, the Berry's phase is theoretically ex- 
tracted from the spinor eigenstate, which are n for single- 
and 2tt for double-layer graphene 3 . From recent stud- 
ies, the tvk Berry's phase in the commonly used form of 
the spinor states is, in fact, related to the pseudospin 
winding number n of a particle as it travels in a loop 
in /c-space which encloses the Dirac point [35J. These 
values have been measured through magneto-transport 
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experiments [4] E] that are typically neither capable of 
measuring Berry's phase of a specific electron bandstruc- 
ture nor free from spin/ valley degeneracy of the elec- 
tron bandstructure of a system under study. Addition- 
ally, this method requires a strong magnetic field, which 
breaks time-reversal symmetry in graphene. Meanwhile, 
the signs of hopping integrals between non-equivalent or- 
bitals for graphene (graphite) have only been determined 
by ab initio calculations J5], e. g. , using maximally local- 
ized Wannier functions 10 . Since the sign of hopping 
integral depends on the characteristics of the localized 
orbitals and the interaction between them, it is crucial in 
determining the electron bandstructure within a tight- 
binding formalism. However, the absence of methodol- 
ogy has led to use the signs following the well-known 
Slonczewski-Weiss-McClure model [TTJ [12] without ex- 
perimental verification for the past few decades. Ad- 
ditionally, the sign of hopping integral between non- 
equivalent orbitals has never been determined experimen- 
tally for any material. 

Given the high momentum-resolving power of angle- 
resolved photoemission spectroscopy (ARPES), ARPES 
is an ideal candidate to solve the above issues on the 
determination of quantum phases. For example, the 
phase difference between the matrix elements describ- 
ing two different optical transitions at the (110) sur- 
face of platinum was extracted from a combined study 
of a spin-resolved ARPES experiment and a theoretical 
model [13] . Also, ARPES has been employed to study the 
characteristics of the spinor eigenstates in graphite |14j 
and graphene |15j . which revealed an interference ef- 
fect between photo-excited electrons [El [15]. However, 
these theoretical studies, within a tight-binding formal- 
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ism, have incorrectly treated the interaction Hamilto- 
nian, which is the key part in the photoemission process 
as it describes the interaction between photons and elec- 
trons. Moreover, it has not been clear how Berry's phase 
enters in ARPES intensity and the sign of hopping in- 
tegral has only been speculated without any comparison 
with experiments [15] , which naturally leads to incorrect 
values. 

Here we report that ARPES can indeed provide in- 
formation on these quantum phases, e. g. , the Berry's 
phase and the sign of hopping integral between non- 
equivalent orbitals. The phase factor in the spinor eigen- 
state of graphene [3] gives rise to strong intensity vari- 
ation around a constant energy contour. Upon rotat- 
ing light polarization by 7r/2, we found that graphene 
with a Berry's phase of n -k (n = 1 for single- and 
n = 2 for double-layer graphene) exhibits the rotation 
of ARPES intensity maxima by 7r/n, which gives impor- 
tant advantages compared to the conventional magneto- 
transport method [U [5|. Additionally, we found that 
full polarization-dependence of ARPES signal reveals the 
sign of hopping integrals in both single- and double-layer 
graphene (graphite can also be understood), e. g. , j' a > 
and "f[ > 0, which is the first experimental determination 
of the sign of hopping integral between non-equivalent or- 
bitals for any material by any method. 



graphene layers at each position, the electron reflectiv- 
ity is plotted as a function of electron kinetic energy, as 
shown in Fig. [ljL) , where the number of dips is the same 
as the number of graphene layers. Regions 1,2, and 3 in 
Fig. 1(a) show 1, 2, and 3 dips, respectively, correspond- 
ing to single-, double-, or triple- layer graphene, respec- 
tively. These regions are painted in black, white, and 
grey, respectively, in Fig. [ljc). The fractions of regions 
in the sample covered by different numbers of graphene 
layers were determined from the areal fractions of differ- 
ently colored regions in Fig. [ljc) . In particular, we find 
that the double-layer graphene sample contains ^74 % 
of double-layer and ^22 % of single-layer graphene. 

III. EXPERIMENT 

We have performed polarization-dependent ARPES 
experiments on single- and double-layer graphene at 10 
K using a photon energy of 50 eV at beam-lines 10.0.1 
and 12.0.1 of Advanced Light Source at Lawrence Berke- 
ley National Laboratory. In Figs. 2(a) and 2(b), we show 
the typical geometry of ARPES experiments: a beam of 
monochromatized light with energy hoj and polarization 
vector e is incident on a sample, resulting in the emission 
of photoelectrons in all directions. The polarization vec- 
tor of light is referenced with respect to the sample nor- 



II. SAMPLE PREPARATION 

Single- and double-layer graphene samples were 
grown epitaxially on n-doped 6_ff-SiC(0001) surfaces by 
electron-beam heating, as detailed elsewhere [IB]. An SiC 
sample was mounted in a prep-chamber with a base pres- 
sure of 5xl0 -10 Torr to remove a thick oxide layer from 
the sample by heating at 600 °C for a few hours. The 
clean sample was then transferred to a custom-designed 
chamber equipped with low-energy-electron microscopy 
(LEEM) with a base pressure of 2X10" 11 Torr and heated 
to 1000 °C under Si flux to improve the surface condi- 
tions for graphene growth. The sample temperature was 
raised to 1400 °C or 1600 °C (determined by an opti- 
cal pyrometer) to make single- or double-layer graphene, 
respectively. 

The surface morphology was monitored in situ dur- 
ing the sample growth by LEEM at the National Cen- 
ter for Electron Microscopy at Lawrence Berkeley Na- 
tional Laboratory. The thickness of fabricated graphene 
samples was determined by LEEM measurements per- 
formed at room temperature following the standard pro- 
cedure [TTJ, [TH]. In particular, the electron reflectivity 
versus kinetic energy curve varies significantly with the 
number of graphene layers providing position-dependent 
measurements on the number of graphene layers. A typ- 
ical bright field image for double layer graphene is shown 
in Fig. [lja) over 4 /imx4 /im range, recorded with elec- 
tron beam of kinetic energy 3.5 eV denoted as the dashed 
line in Fig. ITTb) . In order to determine the number of 




FIG. 1: (a) A LEEM image using an electron energy of 3.5 eV 
over 4/imx4/im range, (b) Reflectivity spectra for the three 
regions (1, 2, and 3) specified in (a), (c) Post-processed image 
of (a) showing the regions covered by single-, double-, and 
triple-layer graphene. (d) A histogram showing the fractions 
of single-, double-, and triple-layer graphene in our sample 
used for double-layer graphene measurements. 
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FIG. 2: (a) X-polarization geometry, (b) Y-polarization geometry. A beam of monochromatic lights with energy Tiuj =50 eV 
and polarization vector e is incident on a sample. The light polarizations in X- and Y-polarization geometries are almost 
parallel to the x and y axes, respectively, (c, d) Measured intensity maps of single-layer graphene at energy E = Ep with 
X- and Y-polarized lights, respectively. Intensity maxima are denoted by white arrows and the electronic band structure of 
single-layer graphene is drawn in the cartoon, (e) Constant-energy ARPES intensity maps for single-layer graphene at Ef with 
X- and Y-polarized light, (f) The angle-dependent intensity profiles of single-layer graphene are obtained by integrating the 
constant-energy intensity map along the radial direction from the Dirac point, in which the solid and dashed lines denote the 
experimental data for X- and Y-polarized lights, respectively. The angle 6 is measured from the +k x direction. The plotted 
quantities are with respect to the intensity minimum. 



mal. In the experiment presented here, two different ge- 
ometries were employed as shown in Figs.[2ja) and[2jb). 
In one geometry shown in Fig. [2ja) , the polarization of 
light is almost parallel to the x axis, while in the other 
shown in Fig. [2jb) to the y axis; hence, we define these 
two geometries as X- and Y-polarization, respectively. 
These geometries have the advantage with respect to the 
conventional s- and p-polarizations used in previous stud- 
ies [HI [5D] , to measure the whole two-dimensional vari- 
ation of the intensity maps around a singular (Dirac) 
point and not just the intensity distributions along two 
characteristic lines in momentum space. This aspect be- 
comes particularly relevant for some experimental condi- 
tions, e. g. , photon energy and sample orientation (i. e. , 
the mixture of light polarizations), when the intensity 
maps (or initial electronic states) are neither symmetric 
nor anti-symmetric with respect to the reflection plane. 
Under this condition in fact, the conventional notations 
would not give appropriate information on the symmetry 
of the initial states. 

Figures [2^c) and [2jd) show the experimental photo- 
electron intensity maps at the Fermi level, E-p, versus the 
two-dimensional wavevector k for single-layer graphene, 
for the two polarization geometries. Here, Ep is 0.4 eV 
above the Dirac point energy En pT] EH [22]. The 
main feature in the intensity maps of both geometries 
is an almost circular Fermi surface centered at the K 



point as shown in Figs. W(c) andlSJd), as expected for 
a conical dispersion. This is in good agreement with a 
recent polarization-dependent ARPES study on single- 
layer graphene when using photons with energy lower 
than 52 eV E3] . Surprisingly we find that the the angular 
intensity distribution is quite different for the two polar- 
izations: for the X-polarization geometry, the minimum 
intensity position is in the first Brillouin zone, whereas 
for the Y-polarization geometry, the maximum intensity 
position is in the first Brillouin zone, suggesting a it rota- 
tion of the maximum intensity in the k x -k y plane around 
the K point upon rotating the light polarization by 7r/2, 
from X to Y (see the white arrows). 

The exact rotation angle is extracted from the di- 
rect comparison between the raw momentum distribution 
curves (MDCs) in Fig. [2je) and the angular dependence 
of the photoelectron intensity maps integrated over the 
radial direction around the K point in Fig. WU) ■ There 
is an overall shift of the intensity maxima (minima) by 
~ 7r upon changing the light polarization from X (black 
solid line) to Y (black dashed line), although the latter 
appears to be slightly shifted by ~ 7r/10 with respect to 
7r. As we will show later, this is due to the presence of 
a finite polarization component along the k x direction in 
our experimental geometry. 

We note that not only the angular position of the 
intensity maximum, but also the absolute value of it 
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FIG. 3: (a, b) Measured intensity maps of double-layer graphene at energies E — Ed +0.25 eV (= Ep) and E = En — 0.95 eV. 
Intensity maxima are denoted by white arrows and the electronic band structure of double-layer graphene is drawn in the 
cartoon, (c) Measured intensity maps of double- layer graphene at energy E = Ep. (d) The angle-dependent intensity profiles 
of double-layer graphene at energy E = Ep, in which the solid and dashed lines denote the experimental data for X- and 
Y-polarized lights, respectively. The plotted quantities are with respect to the intensity minimum. 



changes upon changing light polarization. The maxi- 
mum intensity ratio from experiments is X-polarized/Y- 
polarized=21.4. However, this number itself is not very 
meaningful, because the measured ARPES intensity is 
affected by the difference in the experimental geometries 
for X- and Y-polarized lights (the difference in the out- 
of-plane component of light polarization, photon flux per 
unit area, etc., which are the factors that cannot be 
controlled to be the same in different experimental ge- 
ometries). On the other hand, the ratio from our the- 
ory that will be discussed later provides X-polarized/Y- 
polarized=0.83, assuming that the experimental param- 
eters for two geometries are the same except for the in- 
plane light polarization. 



A similar study on Bernal stacked double-layer 
graphene reveals a strong and complicated momentum- 
, band-, and polarization-dependence as shown in 
Figs, ^a.) and[3jb), that is qualitatively different from 
that of single-layer graphene. Like single-layer graphene, 
the double-layer sample is slightly n doped [25], there- 
fore, only three of the four n bands are occupied and 
hence detectable with ARPES as shown in the cartoon of 
Fig. [3]ja). The most prominent feature is that, when the 
light polarization is changed from X to Y, the maximum 
intensity positions around the K point in the k x -k y plane 
are rotated by ~ ir/2 (see white arrows at Ep). This is 
in striking contrast with the single-layer case where the 
rotation is ~ n as seen from the raw MDCs in Fig. gc) 
and the photoelectron intensity maps integrated over the 
radial direction around the K point shown in Fig. [3jd). 
Due to trigonal warping effects [24], however, the rota- 
tion for higher-energy states is not exactly 7r/2 as shown 
in Figs.[3ja) and|3^b). 



IV. THEORETICAL ANALYSIS 



To the best of our knowledge, the only models in the 
literature describing the polarization-dependence of the 
ARPES intensity in graphite [H] and graphene [TS] are 
substantially different from our results. Previous studies, 
in fact, predict a small polarization-dependence |14j and 
no polarization-dependence |15j of the photoelectron in- 
tensity maps, respectively. Therefore, to be able to repro- 
duce our experimental findings and understand what lies 
behind this nontrivial polarization dependence, we have 
developed a new model. In particular, we first consider 
the Hamiltonian using the tight-binding model based on 
the p z orbital of each carbon atom using two parame- 
ters: t and ti for the in-plane nearest-neighbor (A-B or 
A'-B') and the inter-layer vertical (B-A') hopping inte- 
grals, respectively, as schematically drawn in Fig. |4|a). 
The parameters to and t\ correspond to — j' Q and 7]7 re- 
spectively, in the well-known Slonczewski- Weiss- McClure 
(SWMc) model [TT][T2]. In our calculation, we have used 
|*o | = 3.16 eV and |ti| = 0.39 eV, which are the values in 
Table II of Griineis et al. [9] , but we do not fix the signs of 
them. Note that all four possible choices of the signs give 
exactly the same electron energy band structure within 
this two-parameter tight-binding model. 



With this setup, the tight-binding Hamiltonian of a 
double-layer graphene for two-dimensional wave vector 
k = (fc XI k y ) using a basis set composed of Bloch sums of 
localized orbitals on each of the four sublattices (A, B, 
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(b) 



FIG. 4: (a) Schematic of single- and double-layer graphene. 
(b) The Brillouin zone. Here, bi = 6(0,1), b2 = 

are the three vectors 



and ba = b 



(*-») 



connecting the in-plane nearest neighbor atoms where the 
inter-carbon distance b = 1.42 A, and the lattice constant 
is a — \> r 3b. The positions of the K and K' points are (|^ , 0) 
and (— fj,0), respectively. 



A', and B') reads 

/ t g(k) \ 

t a g*(k) t l e-' lk * d 

t 1 e ik - d t g{k) 

V t Q g*(k) / 



-^double s ) 



Here. 



(k) = exp(ik • bi 



(1) 
(2) 



with bj's defined as in Pig. Ha), and 

^^ e * R A0(r-R A ), (3) 



,ikR B/ 



(4) 



etc. We note that often the k dependence of the ba- 
sis function is suppressed in the spinor notation for sim- 
plicity. In Eq. ([I]), we have considered a phase differ- 
ence e fez d arising from the finite inter-layer distance d 
and the perpendicular component of electron wave vec- 
tor k z . Here, k z is not part of the Bloch momentum, 
but the z component of the photoelectron wave vector, 
which is determined by the photon energy. This quantity 
plays a crucial role in determining the ARPES intensity 
of double-layer graphene as will be discussed later and 
also of multi-layer graphene as previously reported |25j . 

The additional interaction Hamiltonian coupling to 
electromagnetic waves of wavevector Q for a double-layer 
graphene H^ uhlc can be obtained by using the velocity 

operator v — r,H® ouhlc /ih, where f = ih (V k , dk z ) in 
the k-representation and h is the Planck's constant, as 
— | A • v [25] [e is the charge of an electron, c is the speed 
of light, and the external vector potential is given by 
A(r,t) = A Q e'( Q - r -^)] [28] , i.e., 



^doub le (k,Q) = -^ c A Q - 



/ i V k3 (k) \ 

i V k5 *(k) -idt ie -' lk ° d z 

idhe lk - d z ioV k s((k) 

V t V k .g*(k) / 



(5) 



The transition matrix elements in Eq. ^ are those taken 
between basis functions of Bloch sums of p z orbitals 
of wavevectors k + Q and k. Equation ^ is valid 
when 1/|Q| is much larger than the distance between 
the nearest-neighbor atoms b, i. e. , when the variation 
in A(r, t) over a length scale of b is much smaller than 
A(r, t) itself. We shall eventually take the Q — > limit 
in our discussion because the momentum of light is neg- 
ligible for photon energies considered. 

For single-layer graphene, performing a similar type of 



analysis, we obtain 

-^single s ) 

and 



ifi£ gle (k,Q) = -~A Q 



t g(k) 
t g*(k) 



(6) 



* V k5 (k)\ 
i V k5 *(k) 



(7) 



G 



H lnt is critical to explain the polarization depen- 
dence of Ik, because it describes the interaction be- 
tween electrons and photons. The lack of polarization 
dependence in previous studies [TH [T5] is indeed due 
to the way in which H lnt is incorrectly treated. In 
one case [15] . the light interaction is completely ne- 
glected by setting H mt = 1, while in the earlier study on 
graphite [2], the velocity operator v is replaced by the 
momentum p/mo = — i/iV/mo, where h is the Planck's 
constant and mo the free-electron mass. This replace- 
ment works [271 HE] only when the Hamiltonian is lo- 
cal, whereas a tight-binding Hamiltonian, which has been 
used in the previous studies [HI [T5] as well as our study, 
is intrinsically non local. The experimental finding of 
a strong polarization dependence of 1^ in Figs. [2] and [3] 
clearly shows the need for a more complete theoretical 
treatment. We have developed a theory using the widely 
adopted tight binding model with one p z -\ike localized 
orbital per carbon atom, but employing the appropriate 
interaction Hamiltonian with the velocity operator. A 
very good agreement between our model and the exper- 
imental results is obtained for all polarizations and for 
both single- and double-layer graphene, when we com- 
pare Fig. [2] with Fig. [5] and Fig. [3] with Fig. [6] as will be 
discussed later. 

In order to understand what lies behind the observed 
non-trivial and unexpected wavevector-dependent pho- 
toelectron intensity Ik, we need to calculate the abso- 
lute square of the transition matrix element Af s k = 
(/k+cjl H int (k,Q) \ip sis ), where \i/j sk ) is a single- or 
double- layer graphene eigenstate with s = ±1 the band 
index, l/k+cj) is the plane-wave final state projected onto 
the p z orbitals of graphene [both IV'sk) and l/k+cj) are 
expressed using the basis set of Bloch sums of local- 
ized p z orbitals at sublattices A and B in Fig. [5]ja)] and 
H' nt — — I A • v [26], which should not be neglected in 
photoemission process [T5J. The use of a projection of 
the final plane-wave state onto the Bloch sum, which - 
when using plane-waves basis - is effectively composed 
of multiple plane- waves [29] . allow to explain the non- 
trivial polarization dependence of the ARPES intensity 
distribution in Figs. 2(c) and 2(d). Since the polariza- 
tion of A is in the x-y plane, the projection of |/k) onto 
the cr-states of graphene will result in zero contribution to 
the transition matrix elements and hence are neglected in 
this analysis. For simplicity of notation, and without any 
loss of generality, in the rest of this section, we shall take 
the limit of Q ->• and denote iT nt (k) = iJ int (k, Q) and 
l/k) = l/k+cj)- For single-layer graphene, we may use 



and for double-layer graphene 



I/O = 



(8) 




(9) 



For photons (with energy « 50 eV) used in the exper- 
iment, fc z of the planewave final state is much larger 
than the variation of two-dimensional Bloch wavevector k 
around a single Dirac point, leading to only a small vari- 
ation in fc z with any change in k around a single Dirac 
point. For light with a nonzero polarization component 
along the z direction, it will only give rise to an addi- 
tive isotropic term to the photoelectron intensity that is 
independent of the in-plane polarization of the light. 

Now, let us consider the case where k is very close 
to the Dirac point K as shown in Fig. |4|b), and define 
q = k K (|q| < |K|). According to Eq. 



7(q + K) 



V3 



(10) 



where a is the lattice parameter, 
graphene, therefore, 



For single-layer 



Hl gle (q + K) 



— — at (q x a x + q y a y ) , (11) 



and 



single 



(q + K) 



V3e 
~2hc 



at (A 0x a x + Ao y a y ) , (12) 



where a x and a y are the Pauli matrices. The energy 
eigenvalue and wavefunction of Eq. (11) are given by 

E s k = ^o|i |s|q| and 



1 



g-iSq/2 

-sgn(i ) s e ie i/ 2 



(13) 



respectively, when q is the angle between q and the +k x 
direction. Using Eqs. dsb, (12 1, and (13), the transition 
matrix element is given for light polarized along the x 
direction by 



AI 



x — pol 



exp(-z0 q /2) " sg n (*o) s exp^q/2) • (14) 



It follows that for s — +1 (states above the Dirac point 
energy), 



\M 



IE-poll 2 

+ lk I 



and 



\M 



a:— pol i2 
+ lk I 



sm 2 (V2) 



cos 2 ((9 q /2) 



(15) 



(16) 



with tg > and to < 0, respectively. 

Similarly, for light polarized along the y direction, the 
transition matrix element is given by 



Mf" po1 ~ expH V 2 ) + s gn(*o) a exp(i0 q /2) 



(17) 



It follows that for s 
energy), 



T (states above the Dirac point 



\M 



y-poli2 
+ lk I 



! (V2) 



(18) 
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FIG. 5: (a, b) Calculated intensity maps of single-layer graphene for X- and Y-polarized lights, respectively. The insets are 
the results of calculations [TJ] using the simplified momentum operator instead of the correct velocity operator. An arbitrary 
energy broadening of 0.10 eV has been used. Intensity maxima are denoted by white arrows, (c, d) The angle-dependent 
intensity map of single-layer graphene for X- and Y-polarized lights, respectively, in which the solid black, solid red, and dashed 
blue lines denote the experimental data, the calculated results obtained by assuming the actual light polarization used in the 
experiment, and the calculated results obtained by assuming perfectly Y-polarized light, respectively. The theory results shown 
in (a) and (b) have adopted the light polarization used in the actual experiment shown in Figs. |ja) and^b)). The plotted 
quantities are with respect to the intensity minimum. 



and 



im^ o1 i 2 



sin 



V2) 



(19) 



with to > and to < 0, respectively. 

In both cases, we can explain the rotation of the inten- 
sity maximum in the photoelectron intensity map around 
the K point by it upon the change from X- to Y-polarized 
light. Comparing Eqs. (15 1 and (16) with Eqs. (18) 
and (|19), irrespective of the sign of t , the maxima of 



the photoemission intensity map of single-layer graphene 
is rotated by 7t when the light polarization is rotated 
by 7t/2, in agreement with experiment shown in Fig. [2] 
Moreover, the theoretical results with to < shown in 
Figs.[5ja) and[5](b) agree with the measured angular spec- 
tral intensity shown in Figs. f2^c) and [5Jd); especially, 
the choice of \t \ = 3.16 eV (fitted to previous experi- 
ments [9J) reproduces quite well the salient features in 
the experimentally measured intensity maps. 

This is even more clear from the angular dependence 
of theoretical photoelectron intensity drawn with the red 
solid lines compared to experimental results drawn with 
the black solid lines in Figs. [5jc) and [5](d) . Note that 
experimental intensity maximum for Y-polarized light 
shows additional shift by ~ n/10 in Fig. ffl(f). This ad- 
ditional shift is well understood by a finite polarization 
component along the k x direction. When we assume the 
actual light polarization used in the experiment shown in 
Fig. [2|b) , the theoretical intensity exactly matches with 
the experimental result. On the other hand, when we 
assume the ideal Y-polarization, the intensity maximum 
appears at it as shown with the blue dashed line in in 
Fig.[5|d). Therefore, we determine through experiment 
that the inter-orbital hopping matrix element to between 
two in-plane nearest-neighbor carbon p z orbitals is nega- 
tive; we will come back to this point later. 

Recent theoretical study on the matrix element in 
single-layer graphene |23j has found that, in order to de- 



scribe the matrix element for Y-polarized light from first- 
principles calculations using plane- wave basis, one needs 
multiple plane-wave components for the final photo- 
emitted electron state. Since we consider a projection of 
the final state onto the Bloch sum, which - when using 
the plane-waves basis - is effectively composed of multi- 
ple plane- waves [3pJ, our approach can explain the non- 
trivial polarization dependence of the ARPES intensity 
distribution in Figs. 2(c) and 2(d). Additionally, our re- 
sult obtained by using photons with energy 50 eV is in 
good agreement with the recent study [33] based on first- 
principles calculations using photons with energy lower 
than 52 eV. This suggests that, at photon energy below 
52 eV and only when the correct interaction Hamilto- 
nian is employed, the projection of the final state onto 
the tight-binding Bloch sum may describe the true final 
state qualitatively. However, since our theoretical frame- 
work is within tight-binding formulation, and not from 
first-principles calculations with plane- wave basis, con- 
vergence tests with respect to the number of plane-waves 
and the character of the true final state are beyond the 
scope of this paper and, in fact, has been done in a recent 
study [25] . 

In the case of double-layer graphene, for simplicity 
of the analysis, we confine our discussion to the inner 
parabolic bands (bands 1 and 2 in the cartoon in Fig.|3|, 
although we considered all the four bands in our the- 
oretical calculations shown in Figs. |6][a) and [6jb) . In 
double-layer graphene, for electronic states with energy 
\E\ -C of the Hamiltonian in Eq. (IT]), the energy 
and wavefunction are given by -E s k ~ sfrrq 2 /2m* where 
m* = h 2 \t 1 \/2tl and 



1 

71 



y -sgn(ii) se 



(20) 




FIG. 6: (a, b) Calculated intensity maps of double-layer graphene for X- and Y-polarized lights, respectively. An arbitrary 
energy broadening of 0.10, 0.10, and 0.15 eV have been used for bands 1, 2, and 3, respectively. Panels denoted by 'Ideal DC 
show results obtained by considering a double-layer graphene alone and those denoted by 'SG+DG' show results considering 
the contribution from some single-layer fraction of the sample as well (see text). Intensity maxima are denoted by white 
arrows, (c, d) The angle-dependent intensity map of single-layer graphene for X- and Y-polarized lights, respectively, in which 
the solid black, solid red, and dashed green lines denote the experimental data, the calculated results obtained by considering 
the contribution from single-layer graphene portion of the sample, and the calculated results for ideal double-layer graphene, 
respectively. The theory results shown in (a) and (b) have adopted the light polarization used in the actual experiment as 
shown in Figs. [21 a) and[2[b). The plotted quantities are with respect to the intensity minimum. 



respectively [3D]. The phase difference e lkld between the 
two graphene layers in Eq. ([!]) appears here as well. Using 
Eqs. (J5j), ([9]), and (20), the transition matrix element is 
given for light polarized along the x direction by 



M„f- pol ~exp(- 



It follows that for s = 
energy), 



- sgn(fi) s exp(i# q + ik z d) . (21) 
-1 (states above the Dirac point 



I AT 



x— pol|2 
-lk I 



and 



x— pol 1 2 
+ lk I 



2 (0 q + fc z d/2) 



! (# q + M/2) 



(22) 



(23) 



with t\ > and t\ < 0, respectively. The perpendicular 
component of the wavevector k z reads |31j . 



fc z = V 2r n e {E KE + V innel )/H 2 - k 2 



(24) 



where -Eke is the kinetic energy of the photo-electron 
and Mnner the inner potential. Note here that k is a 
two-dimensional Bloch wavevector, i, e. , it has no out-of- 
plane component. The inner potential has been measured 
for graphite by analyzing the energy dispersion at normal 
emission (i. e. , k = 0) [31] , According to Eq. (24), k z d rj 
3.9tt. 

Similarly, for light polarized in the y direction, the 
transition matrix element for a double-layer graphene is 



given by 



M 



y-pol 



exp(~i# q ) + sgn(ii) s exp(z6* q + ik z d) . (25) 



It follows that for s = +1 (states above the Dirac point 
energy), 



\M 



y-poli2 
lk I 



cos 2 (6» q + A; z £i/2) 



and 



sin' 



\e^ + k z d/2) 



(26) 



(27) 



with t\ > and t\ < 0, respectively. 

In both cases, we can explain the rotation of the in- 
tensity maximum in the photoelectron intensity map 
around the K point by 7r/2 upon the change from X- 
to Y-polarized light. Comparing Eqs. (22) and (23) with 
Eqs. ([26| and (27), irrespective of the sign of t\, the max- 



ima of the photoemission intensity map of a double-layer 
graphene is rotated by 7r/2 when the light polarization 
is rotated by tt/2, in agreement with experiment shown 
in Fig. [3] If we assume that t\ > 0, which is qualita- 
tively in agreement with experiment shown in Figs. [3^ a) 



and 3[b), I 



x— pol 
+ lk 



1-1 k 



\ M TiT\ 2 - sin 2 (0 q + k z d/2) and 
cx |M^p o1 | 2 - cos 2 (6» q + k z d/2) for the upper 
band (s = +1). Therefore, we have determined that the 
vertical inter-layer hopping matrix element t\ between 
two carbon p z orbitals sitting on top of each other is pos- 
itive; we will come back to this point later. 
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Although this model can overall account for the ex- 
perimental data of double-layer graphene, there is a 
discrepancy: the measured photoemission intensity at 
E = E-q + 0.25 eV along the +k x direction for the X- 
polarization geometry is finite, whereas the theory pre- 
dicts this value to vanish as shown in the £b + 0.25 eV 
map of the Ideal DG in Fig.[6][a). We believe that this dis- 
crepancy arises from the finite size of the light beam spot 
(~80x40 A*m 2 ) which covers not only the double-layer 
graphene portion but also some single-layer graphene, 
as discussed in Fig. [T] In fact, double-layer graphene 
samples inevitably contain a finite amount of single-layer 
graphene fTTJ [TS]. The fraction of single-layer cover- 
age can be obtained by LEEM measurements jXTl IT8] : 
from this analysis shown in Fig. [T] we find that the 
double-layer graphene sample used here contains ~74% 
double-layer and ^22% single-layer graphene. When the 
theoretical photoelectron intensity maps of single- and 
double-layer graphene are correspondingly weighted and 
averaged, the results denoted by SG+DG in Figs. [6^a) 
and |6|;b) are in excellent agreement with the experimen- 
tal data shown in Figs.[3]Ja) and[3]Jb). This is even more 
clear from the angular dependence of the photoelectron 
intensity maps shown in the red solid lines in Fig. |6jc) 
and[6^d). Note that the presence of single-layer graphene 
does not affect the results for the Y-polarization, which is 
obvious when we compare the red solid and green dashed 
lines in Fig. 6](d) , because both the intensity maxima of 
single- and double-layer graphene occur near 9 = tt. 



of light polarization by 7r/2, i. e. , the ir/n rotation of Ik- 
This prediction is exactly realized in our experimental 
results. 

The power of our method is that it can be extended 
in a straightforward way to other materials with Berry's 
phase j3 — nn (not necessarily tt or 2ir). In this case, 
the photoemission intensities for X- and Y-polarization 
geometries are given by /^~p o1 qt sm 2 ( n # q /2 + ip) 
and /£~ po1 cx cos 2 (n # q /2 + ip), where tp is a system- 
dependent constant. The important feature is that the 
rotation of light polarization by 7r/2 results in a rotation 
of intensity maxima by 7r/n for a system with f3 = rnr re- 
gardless of the constant (p. Thus, we have demonstrated 
here that the Berry's phase can directly be measured 
from polarization-dependent ARPES. 

Unlike methods based on magneto-transport experi- 
ments [H [5] , our new method has three important advan- 
tages, (i) The Berry's phase of a specific electronic band 
can be measured, because ARPES has the angle-resolving 
power and also because one can set up a tight-binding 
Hamiltonian focussing on only the electronic states of in- 
terest: those two results can directly be compared with 
each other, (ii) Due to the angle-resolving power, the 
measured result is free from valley degeneracy for the 
case of graphene. (iii) Our method docs not need electric 
gating, which is essential for the transport measurements. 



B. The sign of inter-orbital hopping integrals 



V. DISCUSSION 
A. Berry's phase 

We have shown that, when the light polarization is 
rotated by 7r/2, the maximum intensity position in 1^ 
in the k x -k y plane of single- and double-layer graphene 
is rotated by 7r/n, where n — 1 and n = 2 for single- 
and double-layer graphene, respectively. The physical 
meaning of these rotations, whose origins rests on the 
phase factor exp (±i n 6* q /2) of the sublattice amplitude 
of the wavefunctions [24] , becomes clear upon a complete 
circulation of q around the Dirac point K (# q — >• 6> q + 
271"), which directly gives a Berry's phase /3 = n A# q /2 = 
nix [24]. Recently the Berry's phase interpretation of 
n [3H5] has been given by a different interpretation in 
terms of pseudospin winding number [35] . 

The fact that n# q /2 enters in 1^ in the form 
of either sin 2 (n # q /2 + ip) or cos 2 (n# q /2 + (p) with 
some constant tp, demonstrates that the matrix el- 
ements directly contain information on the Berry's 
phase. The rotation of light polarization gives 
an additional phase ir/2 to the phase factor, i.e., 
exp (±i n q /2) -> exp (±in6» q /2 + ztt/2). Thus, 
the photoelectron intensity /k is modified accord- 
ingly from sin 2 (n # q /2 + tp) — > cos 2 (n# q /2 + tp) or 
cos 2 (n# q /2 + <p) —> sin 2 (n6* q /2 + ip) upon the rotation 



Another important finding of our study is that we can 
directly extract, for the first time, the sign and the ab- 
solute magnitude of the inter-orbital hopping integrals 
(IOHIs) between non-equivalent localized orbitals of a 
tight-binding Hamiltonian from experiment. Until now, 
in fact, the sign determination of an IOHI has resorted 
no to any experimental method, but to ab initio calcu- 
lations, e. g. , using maximally localized Wannier func- 
tions [10] . In order to understand an ambiguity related 
with the experimental sign determination, we take the 
simplest one-dimensional example, and extend the dis- 
cussion to a more complicated tight-binding model of 
graphitic systems than the one described previously. We 
consider simple tight-binding models having s-like local- 
ized states, the values of which are all positive in real 
space (we can arbitrarily set this gauge without losing 
generality.) If there is only one localized orbital per unit 
cell in a one-dimensional tight-binding model as drawn 
in Fig. [7ja), the energy bandstructure varies with the 
sign of the IOHI t' as shown in Fig. [7|b). Hence, the 
sign of IOHIs between "equivalent" orbitals can always 
be trivially determined [33] . 

However, if we consider the case where there are two 
non-equivalent localized orbitals <f> s and <fi' s whose value 
in real space is positive and if we denote the IOHI be- 
tween the nearest neighboring orbitals by t" as drawn in 
Fig. [7^c) , the actual band structure is invariant when we 
change the sign of t" as shown in Fig. [7^cl) . Therefore, 
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FIG. 7: (a) Schematic of a one-dimensional crystal having 
s-type orbital per unit cell. The nearest neighbor hopping 
integral is i! . (b) Calculated electron energy band structure 
of the system depicted in (a) with different choices for the 
sign of t'. (c) Schematic of a one-dimensional crystal having 
two nonequivalent s-type orbitals per unit cell (i. e. , having 
different on-site energies). We assume that all the distances 
between the centers of the nearest neighbor orbitals are the 
same and hence the corresponding hopping integrals, denoted 
by t" . (d) Calculated electronic band structure of the system 
depicted in c with different choices for the sign of t". 



even when the actual electronic band structure is empiri- 
cally determined, the sign of t" cannot be determined. In 
general, an empirical tight-binding model with more than 
one orbital per unit cell has this sign ambiguity problem 
for IOHIs between "non-equivalent" orbitals, thus pre- 
venting an experimental measurement of IOHI from just 
the energy band dispersions. 

We could understand this degeneracy as follows. If wc 
denote the Bloch sums of the two localized orbitals by 
<pi (k) and 02 (k), the tight-binding Hamiltonian H using 
this basis set reads 



H(k) 



h(k)\H\MV) 
ia(k)|ff|0i(k)) 



(<Mk)|ir|<Mk)) 



(28) 



Note here that the matrix elements (0j(k)| H |<^(k)), in 
which i.j £ {1,2}, involve not only the on-site or nearest- 
neighbor hopping processes but also all the other possible 
hopping processes. Now, it is obvious that the following 
Hamiltonian H'(k) has exactly the same eigenvalues as 
H(k): 



H'(k) = 



6 1 (k)| J rT|0 1 (k)) 
6 2 (k)|iJ|<Mk)) 



b 2 (k)\H\MV) 



(29) 

What we have done by going from H(k) to H'(k) is 
to change the sign of the IOHI between the two non- 
equivalent localized orbitals. In fact, the two matri- 
ces -ff(k) and H'(k) are related by a unitary trans- 
form, which does not change the eigenvalues of a matrix, 
H'(k) = U^H(k)U with 



U = 



1 

-1 



(30) 



On the other hand, if we change the signs of the diagonal 
matrix elements (0i(k)| H |0,(k)), in which I e {1,2}, 
we get a different eigenvalue spectrum. Thus, there is 
no ambiguity in the sign of the IOHI between equivalent 
orbitals. This simple example illustrates that one cannot 
determine the signs of the IOHIs between non- equivalent 
orbitals just by looking at the measured electron energy 
bandstructure. 

The tight-binding model for double-layer graphene 
that we used in our calculations is based on the p z or- 
bitals of carbon atoms with two parameters: to an d ti for 
the nearest-neighbor in-plane and the vertical inter-layer 
hopping integrals, respectively. The parameters to and t\ 
correspond to — 7 and 7^ in the well-known Slonczewski- 
Weiss-McClure model, respectively [TTJ [T^] , and we have 
used |t | = 3.16 eV and |ti| = 0.39 eV 0. The photoe- 
mission intensity map in the k x -k y plane is strongly de- 
pendent on the signs of both t and ti as shown in Fig. [8] 
Because the four different choices of the signs produce ex- 
actly the same electronic band structure, it has not been 
known from experiments which choice of the signs of the 
IOHIs is physically correct, although the absolute values 
have been experimentally estimated [MJ [35] . 

This sign-ambiguity problem still exists even when 
we include more complicated hopping processes in the 
model, especially the non- vertical inter-layer hopping in- 
tegrals 73 and 74, there still exist unitary transforms 
that leave the energy eigenvalues unchanged. In a tight- 
binding model having four hopping integrals (7q, 7J, 7 3 , 
and 7 4 ), the following four different sets of parameters 
give exactly the same electron energy bandstructure: (7q, 

Txj 7s; 74)> (7o! -7x5 -7 3 5 -74)5 ("7o5 7x5 7 3 > -"fi)> and 
(— 7q, — -f[, —73, 7 4 ), assuming that the first set is com- 
posed of the values currently accepted and used when the 
SWMc model is considered. In principle, there should 
be eight different sets of parameters giving the same en- 
ergy bandstructure; however, from our knowledge that 
the nearest-neighbor intralayer hopping integrals in dif- 
ferent graphene layers are the same, we have reduced the 
number of candidates to four. This is also the reason why 
we considered only the four cases in Fig. [8] 

We have shown that the choice of to < and t-y > 0, 
i.e., 7 > and j[ > in Fig. Mc), reproduces the 
experimental results shown in Figs. [3ja) and [3]^b) , thus 
experimentally determined the signs of these IOHIs. The 
signs of the IOHIs used for graphite in the conventional 
model [TTJ [12] is indeed correct. Previous theoretical 
study has also tried to determine the signs [15] . but in- 
appropriate theoretical approaches (as previously men- 
tioned) and the lack of full polarization-dependent ex- 
perimental data have led to incorrect speculations, which 
do not agree with the conventional model [TTJ [T^] as well 
as our results. Our method can generally be used to de- 
termine the sign of the hopping integrals in complicated 
materials such as cuprates as well as in simple materi- 
als such as gallium arsenide and one-dimensional crystals 
having two atoms per unit cell. 
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FIG. 8: (a-d) Calculated intensity maps for four different choices of the signs of nearest-neighbor in-plane and vertical inter-layer 
hopping integrals, to and ti, respectively. Case (c) agrees with the experimental results. 



VI. SUMMARY 



ence properties of a huge variety of materials. 



We have shown that ARPES can be used as a power- 
ful tool to directly measure quantum phases such as the 
Berry's phase of a specific electronic band in graphene 
with advantages compared to the interference type of 
measurements [36ff38] which do not give any information 
on the band-specific Berry's phase, and the sign of the 
hopping integral between non-equivalent orbitals, never 
measured for any material before. The experimental and 
theoretical procedures developed here can be applied in 
studying the electronic, transport, and quantum interfer- 
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